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ABSTRACT 

Photospheric magnetic vector maps from two different instruments are used to model the nonlinear force-free 
coronal magnetic field above an active region. We use vector maps inferred from polarization measurements 
of the Solar Dynamics Observatory /Helioseismic and Magnetic Imager (HMI) and the Solar Optical Telescope 
Spectropolarimeter (SP) aboard Hinode. Besides basing our model calculations on HMI data, we use both, 
SP data of original resolution and scaled down to the resolution of HMI. This allows us to compare the model 
results based on data from diff'erent instruments and to investigate how a binning of high-resolution data eff'ects 
the model outcome. The resulting 3D magnetic fields are compared in terms of magnetic energy content and 
magnetic topology. We find stronger magnetic fields in the SP data, translating into a higher total magnetic 
energy of the SP models. The net Lorentz forces of the HMI and SP lower boundaries verify their force-free 
compatibility. We find substantial diff'erences in the absolute estimates of the magnetic field energy but similar 
relative estimates, e.g., the fraction of excess energy and of the flux shared by distinct areas. The location 
and extension of neighboring connectivity domains difl'ers and the SP model fields tend to be higher and more 
vertical. Hence, conclusions about the magnetic connectivity based on force-free field models are to be drawn 
with caution. We find that the deviations of the model solution when based on the lower-resolution SP data are 
small compared to the diff'erences of the solutions based on data from diff'erent instruments. 
Subject headings: Sun: photosphere — Sun: corona — Sun: surface magnetism — Sun: magnetic topology 



1. INTRODUCTION 

Most recent studies dealing with the magnetic struc- 
ture of the solar corona above active regions use dif- 
ferent force-free model app roaches (see recent review of 
IWiegelmann & Sakuraill2012h and base the modeling on pho- 
tospheric vector magnetic field data from either the Helioseis- 
mic and Magnetic Imager (HMI) on board the Solar Dynam- 
ics Observatory (SDO) or the Spectropolarimeter (SP) of the 
Solar Optical Telescope (SOT) on board the Hinode space- 
craft. Such model approaches are used to compensate the lack 
of routine direct measurements of the coronal magnetic field 
v ector. 

iHao et al.] (120121) used SP data as lower boundary condition 
to an optimization approach to analyze the coronal magnetic 
field associated to a white light flare. The modeling suggested 
that the flare originated from sheared and twisted field lines 
with l ow altitudes bridge d by a set of higher magnetic field 
lines. Ilnoue et aH (1201 2b applied a MHD relaxation method 
based on SP data to investigate the buildup and release of 
magnetic twist and suspected the importance of the relative 
handedness of twist ed field lines and the ambient field. 

ISun et al.l (1201 2b|) employed an optimization method based 
on HMI data to model the temporal evolution of the coronal 
field of an active region over five days. The modeling dis- 
played distinct stages of the build-up and release of magnetic 
energy and analyzed the association to changes in the mag- 
netic field. In a subsequent study, I Sun et al.l (l2012al) used 
the same method for a detailed analysis of the field topol- 
ogy during a series of eruptions observed from HMI. Above 
the apexes of cusp-like loops observed in coronal images, the 
modeling result suggested the presence of a coronal null point. 

Comparisons of the outcome of different force-free model 
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algorithms based on the same lower bou ndary conditions have 
been performed in the past too (e.g., iDeRosa et al.l 120091: 
iGilchrist et al.l 120121) . These studies revealed that order-of- 
magnitude estimates of these models based on large enough 
fields-of-view and high enough spatial resolution of the vec- 
tor magnetic field data can be expected to be reliable. 

The model outcome of the same reconstruction algorithm 
based on d ata from differen t instm ments has been studied 
recently by iThal mann et al.l (120121) . using data from HMI 
and the Vector-S pectroMagnetograph (VSM) of the SOLIS 
project (iKeller et al .1120031) . They found agreements of the re- 
sulting force-free models in form of, e. g., the relative amount 
of energy to be set free during an eruption but also found 
a considerable diff'erence in the absolute model energy esti- 
mates. In particular, they found the estimated energy con- 
tent of the VSM model being about twice of that of the HMI 
model. 

In the present study we use vector magnetic field data of 
HMI and SP as an input to the same force-free model and 
compare the results. This is motivated due to the data of these 
two instruments being widely used recently and in the future 
expected to be frequently used as an input for the modeling 
of the coronal magnetic field. We regard this as important 
in order to test the consistency of the model solutions and at 
the same time to give a feeling about the accuracy of the esti- 
mated physical quantities and magnetic field topology based 
on those models. We also investigate the eff'ect which a bin- 
ning of the SP instrument data to a lower resolution has on the 
model outcome. We, however, do not search for explanations 
of the diff'erences of the HMI and SP inversion products itself, 
considering this as to be out of the scope of this study. 



2. DATA AND ANALYSIS METHODS 
2.1. Data Sources 
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The HMI on board SDO (ISchou et al.ll2012l) obtains filter- 
grams at the photospheric Fei 617.3 nm spectral line. The 
full Stokes vector is retrieved from filtergrams averaged 
over about 12 min and inv erted using the Milne -Eddington 
(ME) inversion algorithm of lBorrero etaH (120101) . The 180°- 
azimuth ambiguity of the trans verse field is resolved using 
a Min i mum Energy Alg orithm (iMetcalJI 19941: iMetcalf et~aD 
120061: iLeka et aLri2009l) and the resulting vector magne- 
tograms have a plate scale of 0.5 arc-second. 

In its fast scan mod e, SP as part of the SOT 
(iTsuneta et al. 2008; Sue matsu et al.l 120081: llchimoto et al.l 
120081: IShim izu et ^ |2008|) on board the Hinode spacecraft 
(iKosugi etal.. .2001) obeys a spatial resolution of ~ 0.32 arc- 
second. It observes the Fe i spectral line couplet at 630.15 nm 
and 630.25 nm. Full Stokes profiles are obtained with a 
spectral sampling of 2.25 10"^ nm and a slit scan time of 
1.6 s. The physical parameters from the full Stokes pro- 
files were obtained using the MERLIN ME inve rsion algo- 
rithm (ISkumanich & Liteslll987l:lLites et al.ll2007h . The 180°- 
azimuth ambiguity is resolved in the same way as for HMI 
data. 

2.2. Event Selection and Data Set 

To perform a study as described above simultaneous ob- 
servations from HMI and SP are required. Both, HMI and 
SP vector magnetic field data are available to analyze the 
magnetic structure of active region 11382 on 2011 Decem- 
ber 22. SP scanned this active region from 04:46 UT to 
05:29 UT (i. e., scanned ~ 3 Mm/min from solar east to west). 
The HMI vector map used for this study was retrieved at 
Tree = 05:00 UT, approximately at half of the scanning time 
of SP. Given average photospheric conditions (5~ lOOmT, 
p-lO""^ kgm"^) the Alfven speed is VA~10kms"^ The 
Alfven travel time and thus the evolution of the photospheric 
magnetic field over a characteristic distance in the present 
study (~ 100 Mm) is ~ 1 h. Therefore, we assume that the 
temporal changes of the photospheric field over ~ 20 min, 
during which SP scanned before and after the time when 
HMI recorded, are negligible. Regarding the longitudinal 
magnetic field we assume this as justified since, e. g., also 
I Wang et aP (l2009h found that the average ratio of the longitu- 
dinal magnetic flux d ensities measured by SP and SoHO/MDI 
(IScherrer et al.ll 19951) was not strongly influenced by the evo- 
lution of the photospheric magnetic field during the scanning 
time of SR 

The magn etic field vectors are tran sformed to Heliographic 
coordinates (iGary & Hagyard|[l990h . We correct the SP data 
set for the eff'ect of diff'erential rotation where we use Tree 
as reference time. Given present computational capabilities 
high-resolution data is sometimes binned to a lower resolu- 
tion in order to allow for near real-time magnetic field model- 
ing and/or computational domains of feasible dimensions (see 
e. g., DeRosa et al. 2009). Thus, in the course of co-alignment 
of the HMI and SP data, we also bin the original-resolution SP 
(SPorig) data to the resolution of HMI which involves a 2D lin- 
ear interpolation. The binned SP data is hereafter denoted as 
SPbin data. 

The field-of-view (FOV) used to study active region 11382 
is mainly determined by the area covered by the SP scan 
(~140x85Mm^ centered around S19W07; see Figure [B- 
Within the HMI data which covers roughly ~ 400 x 250 Mm^, 
we define a window equally sized as the FOV of SP and calcu- 
late the cross-correlation between its vertical field component 
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Fig. 1. — (a) HMI and (b) SPbin data. The gray-scale background reflects 
the vertical magnetic field, B^. Black/white arrows indicate the direction 
of the horizontal magnetic field, originating from positive/negative polarity 
where B^ > 20 mT. The dashed rectangle outlines the quiet-Sun region used 
to calculate the Icr uncertainty of the vertical and horizontal magnetic field. 
The black arrow just below the dashed rectangle indicates the length of an 
arrow representing a horizontal field magnitude of 100 mT. 

and that of the SP data. By shifting the position of the HMI 
window we search for the highest cross-correlation to find the 
corresponding HMI sub-field. 

2.3. Uncertainty Estimation 

The noise level of the HMI data is on the order of 
ImT/lOmT for the longitudinal/transverse magnetic field 
(X. Sun, private communication). The average uncertainty 
for both the longitudinal and transverse field, estimated from 
the SP inversion error maps is ^ 1 mT. Besides seeming rather 
low especially for the transverse field these standard error es- 
timates may not be reliable (B. Lites, private communication). 
Thus, we employ a consistent measure of the uncertainty level 
for the data from the two instruments, as described in the fol- 
lowing. 

When investigating the properties of the HMI and SP data 
in the course of the force-free modeling, we only consider 
(i) pixels with values of the vertical field, B^, and horizon- 
tal field, Bh, above the respective 2cr uncertainty levels (2(Tb^ 
and 2(Tfi/^, respectively) or (ii) pixels with values ofB^ > 2crB^ 
and Bh < ^o-Bh ^ the latter in order not to disregard the strong 
vertical fields in the center of the active region. 

The respective 2(T uncertainty levels are calculated for the 
quiet-Sun area, outlined by the dashed rectangles in Figure [H 
where we find: Icfb^ = 3.3/6.7 mT and 2crB,^ =5.2/12.2 mT for 
the HMI/SPbin data. The uncertainty levels for the SPorig 
data are calculated from a quiet-Sun region, equivalent to 
that outlined in Figure [T] Here we find 2crB^ = 8.8 mT and 
2(Tfi/^ = 12.3 mT, i. e. slightly higher than what was found for 
the SPbin d ata. The latter estimates, in fact, conform with the 
findings of lOrozco Suarez et al.l (l2007h who had estimated the 
measurement uncertainties of SP for the quiet-Sun inter net- 
work field strengths and fluxes as < 15 mT. 

2.4. Magnetic Field Modeling 
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Fig. 2. — Relative occurrence of the absolute (a) vertical field, B^, and (b) horizontal field, Bh. The power of B^ and Bh as a function of length scale, /, and wave 
number, k, is shown in panels (c) and (d), respectively. HMI/SPbin data are represented by black/gray lines. Only pixels with values B^ > 2o-b, and Bh,> '^ctb,^ or 
B^ > 2o-B^ and Bh < ^os,^ are taken into account. 



Photospheric polarization signals originate from atmo- 
spheric ]^Yers_whi£hareknown not to be force-free. For in- 
stance, iMetcalf et alJ (119951) showed, using vector magnetic 
field measurements of an active-region magnetic field, that it 
can be considered as to be force-free above 0.4 Mm above a 
photospheric level. Modeling t he interchang ing dominance of 
plasma and magnetic pressure, iGaryl (1200 lb was able to esti- 
mate the height regime of dominating magnetic fields above 
a sunspot/plage region as -0.8 - 200 Mm. Using high- 
resolution SP vector maps, iTiwaril (120 12h found that umbral 
and inner penumbral parts of sunspots may nearly be force- 
free but that sunspots as a whole at a photospheric level might 
not entirely be so. 

Therefore, the inferred HMI and SP magnetic vector maps 
are not force-free consistent and need to be preprocessed 
to achieve suitable force-free consistent boundary conditions 
dWiegelmann & Inhester 2006) . From the vertical component 
of the preprocessed field vector, a potential field is calculated 
and used as start equilibrium and to prescribe the boundaries 
of the cubic computational domain. The bottom boundary 
is replaced by the preprocessed vector field and the set of 
force-free equati ons for the nonlinear case in C artesian co- 
ordinates solved dWiegelmann & Inhesteill201Ql) . A bound- 
ary layer of ~ 10 Mm is introduced towards the lateral and 
top boundaries where the nonlinear force-free (NLFF) solu- 
tion drops to the prescribed boundary field. For our analy- 
sis we discard this layer and only consider the inner (phys- 
ical) ~ 120 X 60 X 70 Mm^ domain and the according bottom 
boundary field. Since this method involves the relaxation of 
the magnetic field not only inside the computational domain 
but also on its bottom boundary we compute a potential field 
from the relaxed lo wer boundary based o n the Fast-Fourier 
method described by lAlissandrakisI (|1981). 

Hereafter, we refer to the 3D NLFF field model based on 
the HMI vector map as an input to as "HMI model". Simi- 
larly, the "SPbin model" and "SPong model" result from using 
the binned and original-resolution SP vector magnetic field, 
respectively, as input data to our preprocessing and force- free 
reconstruction algorithms. We compare the HMI and SPbin 
model in § I3.H43.3| and summarize the eff'ects of binning on 
the model outcome in § 13.41 



3. RESULTS 
3.1. Magnetic Flux and Force-Freeness 



TABLE 1 
Magnetic field related quantities of the vector data 
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Listed are the total unsigned vertical flux, |0^|, the average horizontal field, Bh 
as well as the net Lorentz force components (\Fx\, \Fy\ and IF^I) normalized 
to the magnetic pressure force, |Fp|. Only pixels with values B^ > ^o-b, and 
Bh,> 2o-B,^ or Bz > 2o-B, and Bh < 2o-b,^ are taken into account. 

After performing the data preparation as described in § 12.21 
we are able to investigate how the vertical and horizontal field 
components of the HMI and SPbin vector maps compare to 
each other and to check the force-free consistency. We o nly 
consider data points where the criteria outlined in § 12.31 are 
fulfilled. 

We find on overall stronger vertical fields in the SPbin than 
the HMI data (especially for 52>100mT; see Figure [2^) 
and stronger horizontal field (except for Bh< 10 mT; see Fig- 
ure [2t)). This can be seen also when comparing the area- 
integrated unsigned vertical flux, |0^|, and the average hori- 
zontal field, Bh'. the HMI data hosts -54% of |0zlsPbin and 
-81% of 5/i SPbin (^^^ Table[T]). HMI and SPbin have a compa- 
rable sensitivity on long scales, /, i. e. at small wave numbers, 
k, especially on scales / > 5 Mm (Figure [2]:,d). With decreas- 
ing scale the amount of detected field increasingly diff'ers: SP 
data show considerably stronger fields on smaller scales. 

A necessary condition for a magnetic field to be force- 
free is that the components of the net Lorentz force are con- 
siderably smaller than a characteristic magnitude of the to- 
tal Lorentz force in case of a non force-free magnetic field. 
The latter can be approximated by the magnetic pressure, Fp, 

on the lower boundary (lLowlll984l) . The ratio |/^/|/|i^p| with 
/ = (x, y, z) in Table [T] shows that this conditions are met only 
to a certain degree. The ratios found here agree with the val- 
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Fig. 3. — Power of the absolute vertical field, B^, and horizontal field, Bh, before and after applying a preprocessing to it (black and gray curves, respectively). 
Panels (a)/(c) show the power of B^ of the HMI/SPbin data while panels (b)/(d) show that of Bh as a function of length scale, /, and wave number, k. Only pixels 
with values B^ > 2(Tb, and Bh,> los,^ or B^ > '^ctb, and Bh < '^ctb,^ are taken into account. 



ues found by, e. g.. lMoonet aH (120021) and lTiwaril (120121) and 
we also find |F^| > [ \Fx\, \Fy\ ]. Preprocessing, however, cer- 
tainly improves the force-freeness as it smooths the vertical 
field and, additionally, alters the horizontal field in order to 
minimize the net force and torque and to gain boundary con- 
ditions compatible with the force-free assumption (showing 
ratios of clearly smaller than unity). The eff'ect of smooth- 
ing can be seen when comparing the spectral power of the 
raw and preprocessed HMI and SPbm data (Figure[3]): the sig- 
nal on shorter scales (i. e., for large k) is reduced. On longer 
scales, the power of the vertical field remains the same (see 
Figure [3t,c) but that of the horizontal field is enhanced (more 
pronounced for the SPbm data; see Figure [3j),d). 

The preprocessing leads to a slight reduction of |0^|: |0^|hmi 
is reduced by ~3% and |0^|sp by ~4%. The preprocessing 
also leads to an enhancement of Bh'. Bh^uui is enhanced by 
-4% and Bh^sp^m ^^ enhanced by -30% (see Table [TJ. The 
preprocessed HMI data hosts ~ 55% of l^^lsp and ~ 64% of 
^/i,sPbin ^f the preprocessed SPbm data. In summary, the pre- 
processing only slightly reduces the diff'erence of the vertical 
unsigned flux between HMI and SPbm (which is ~ 45% before 
as well as after preprocessing) but enhances the difl'erence of 
the average horizontal field, ABh, (from ABh -19% before to 
ABh ~ 36% after preprocessing). These preprocessed HMI 
and SPbin data are used as lower boundary condition for the 
NLFF reconstruction. 

As mentioned above, our NLFF modeling algorithm also 
relaxes the magnetic field on the bottom boundary of a cubic 
computational domain. Thus, besides during the preprocess- 
ing, the magnetic field vector on the lower boundary is altered 
also while iteratively seeking for the force- and divergence- 
free field solution in the volume. It's modification to the pre- 
processed lower boundary data, as listed in Table [T] is that |0^| 
is reduced by < 0.1% and Bh is increased by ~ 10%. Thus, the 
HMI NLFF lower boundary hosts ~ 55% of |0Jsp and ~ 62% 
of 5/i SPbin of the SPbin NLFF lower boundary data, comparable 
to the ratio we found for the preprocessed data. 

3.2. Magnetic Energy 

From the 3D model fields, we can estimate the magnetic en- 
ergy content of the potential field (£^pot) and of the NLFF field 
(total energy; £^niff)- An upper limit for the energy which can 
be released (excess energy) is given by AE = E^i^ - E^ot. We 
can estimate the statistical accuracy of our volume-integrated 
energy estimates by adding diff'erent artificial noise models 
to the HMI and SP magnetograms, consecutive application of 



TABLE 2 
Magnetic energies of the 3D force-free models 
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^pot 



AE 



[Xl025 J] 



AE 
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HMI 

SPbin 



3.45 2.80 0.65 0.19 
7.44 5.80 1.64 0.22 
7.31 5.85 1.46 0.20 



Given are the total, potential and excess magnetic energy of the 3D magnetic 
model fields, listed as £^niff, £^pot and AE, respectively. The ratio AE/E^m 
gives the relative amount of excess energy. The statistical error of these esti- 
mates is ~ 1% for Epot and £^nifF and ~ 10% for AE. 

the preprocessing and extrapolation algorithms and compar- 
ison of the resulting energy values. This yields a statistical 
error of ~ 1% for both E^ot and E^m and ~ 10% for AE. 

The estimated absolute potential energy of the HMI model 
is ~ 48% of that of the SPbin model (see TableO. This is a di- 
rect consequence of the HMI NLFF bottom boundary hosting 
only ~ 55% of the unsigned vertical flux of the SPbin bottom 
boundary (see Table [TJ. A similar trend is found for the abso- 
lute total and excess energy of the HMI model. 

However, the relative excess energy is about 20% of the 
total energy in both models (given by the ratio AE/E^i^ in 
Table [21). An excess energy on the order of 10^-^ - 10^"^ J is as- 
sumed to be sufficient for powering C-class flaring, which was 
actually observed for the active region analyzed here on the 
days before. Simila r values related to C-cla s s flaring activity 
were f ound by, e. g. jRegnier & Pries3(l2007h : lThalmann et al.l 
(l2008h and lately bv lGilchrist et al.l(l2012h . 

3.3. Magnetic Field Topology 

Besides comparing a volume-integrated quantity like the 
magnetic energy, we are also interested how the modeled 
magnetic field configurations compare to each other. To do 
so and to ensure that we are looking at the same topological 
structure we look for regions of strong gradients in the mag- 
netic connectivity. They are thought of being linked to the 
creation of strong current concentrations in the solar corona 
and believed to represent the footprint of quasi- separatrix lay- 
ers (iPriest & Demoulinll995h. We Quantify the magnetic con- 
nectivity following iTitov & Homigl (l2002h and calculate the 
squashing degree 2 at a height of ~ 5 Mm above the NLFF 
lower boundary (Figure |4t,b). The squashing degree quan- 
tifies the eccentricity of an elliptical cross-section of a flux 
tube into which a flux tube of initially circular cross-section 
is transformed. Wherever Q is large, the magnetic field con- 
nectivity changes drastically over short distances. According 
to this pattern, we choose a region of interest (ROI) around 
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Fig. 4. — Squashing degree, Q, at a height of ~ 5 Mm above the NLFF lower boundary as determined from the (a) HMI and (b) SPbin model, shown in the range 
(55Mm<x,_y). Field lines originating from locations where Q>100 within the ROI (rectangular outlined area) and which have both footpoints located on the 
lower boundary of the HMI/SPbin model are shown in (c,e)/(d,f). The magnetic configuration when viewed along the surface normal is shown in (c) and (d). 
The gray-scale background reflects the vertical magnetic field, B., of the NLFF lower boundary. Black and white contours are drawn at [-i-50mT,-i-100mT] and 
[-50mT,-100mT], respectively. The view from solar south (along _y) is shown in (e) and (f). Color-coded is the absolute height of the calculated field lines 
in Mm. 
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Fig. 5. — NLFF lower boundary of the (a) HMI and (b) SPbin model in 
the range (55 Mm < x, j). The vertical magnetic field is shown as gray-scale 
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lower boundary are color-coded according to the vertical magnetic field there. 



a clearly distinguishable pattern of high values of Q (rectan- 
gular outline in Figure |4^,b). Though clearly visible in both 
models, the Q-ridge appears more diffuse in the SPbin model 
and its location differs up to Ax ^ 5 Mm for a given y. 

In total, 92/83 magnetic field Hnes in the HMI/SPbin model 
(i) start from locations where 2 > 100 within the ROI and (ii) 
connect back to the lower boundary (see Figure |4};-f). They 



Fig. 6. — Relationship of the length of all closed model field lines, /fl, to 
their apex height (/zmax) for the (a) HMI and (b) SPbin model. Field lines 
are calculated starting from every pixel location in the range (60 Mm < x, y) 
where 5, > 2cr5, and Bh,> '^o'Bh or where B^ > 2o-b, and Bh < '^o'Bh ^^d the 
color-code reflects the value of the absolute vertical field there. 

qualitatively outline two neighboring connectivity domains 
which connect the negative polarity region to its neighboring 
positive polarity surrounding. The locations where the field 
lines connect back to the lower boundary are shown color- 
coded based on the local vertical field in Figured] We assume 
that the field lines we calculated represent thin flux tubes. For 
each flux tube, we choose its cross section at the location from 
where we started the field line calculation, dJ?l, as the size of 
one pixel (i.e., dJ?l~360^ km^). Furthermore, we assume 
the vertical field there determines the flux of the flux tube. 
Summation over the all considered thin flux tubes gives an 
estimate of the total absolute shared flux for the HMI/SPbin 
model, where we find 10^1 = 1.0/1.3 x lO^^Mx. We find a 
lower value of connected flux in the HMI model though more 
closed magnetic field lines are considered due to our selection 
criterion. However, comparable is the relative amount of con- 
nected flux linked by these field lines: \(p[\ comprises ~ 1% of 
the unsigned vertical NLFF lower boundary flux as listed in 
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Fig. 7. — Relative occurrence of the absolute (a) vertical field B^, and (b) horizontal field, B^. The power of 5^ and Bh as a function of length scale, /, and wave 
number, k, is shown in panels (c) and (d), respectively. SPbin data and original-resolution data (SPorig) are represented by gray and black lines, respectively. Only 
pixels with values B^ > 2(Tb^ and Bh,> '^o-Bh or B^ > Ictb. and Bh < '^o'Bh are taken into account. 



Tablell] 

We also recognize differences in the connectivity pattern of 
the field lines as shown in Figure [S] Numerous SPbin model 
field lines connect to the weak-field regions (B^ < 50 mT) be- 
tween the two major positive polarity patches centered around 
(x ~ 65 Mm, J ~ 35 Mm) or to the weak-field surrounding at 
(90 Mm < X, 10 Mm < y). However, none of the selected field 
lines of the HMI model does so, instead they re-enter the 
NLFF lower boundary at locations of strong vertical fields 
(Bz > 50 mT). From Figure |4]:-f one recognizes that a larger 
number of field lines connects to the positive polarity at 
(x,};)- (120 Mm, 5 Mm) and the highest field lines of the 
SPbin model reach up to greater heights than those in the HMI 
model. 

To investigate if the latter represents a general trend or is 
biased due to our restrictive selection of start locations for 
field line calculation, we consider all field lines starting from 
a pixel location on the NLFF bottom boundary in the range 
(60 Mm < X, y) where B^ > 2crB^ and Bh,> ^oSi^ or B^ > ^ctb^ 
and Bh < ^o-Bh • We then find very similar relationships of 
the length of the calculated field lines, /fl, to the height of 
their apex, /z^ax for both models (see Figure [S]): the rela- 
tion is w ell defined and to a first approx imation linear (see 
also, e. g. ISchrijver & Aschwandenll2002|) . The SPbin model 
field lines follow a steeper distribution, i.e., seem to be on 
overall higher. It appears that the field lines carrying the 
strongest vertical fluxes are neither the shortest ones nor the 
longest ones. Instead, it seems that in both models, field 
lines with a length of 20 < /fl ^ 100 Mm and an apex height 
of 5 < /zmax ^ 35 Mm carry most magnetic flux. 

The longest and highest closed field lines in the HMI model 
are found to be ~ 150 Mm and 50 Mm, respectively. This 
is lower than the values found for the longest/highest closed 
field lines of the SP^n model (~ 170 Mm/60 Mm). It is not 
surprising that the SPbin model field lines tend to be longer and 
reaching higher up in the model atmosphere when looking at 
the magnetic field distribution on the lower boundary. There, 
we find an average of < BJBu > = 0.8/0.9 for the HMI/SPbin 
model, i. e., the SPbin model field lines are on average more 
vertical. 

3.4. Effects ofSP data binning on the modeling 

The binning of the original-resolution SP data (SPorig) to 
the resolution of HMI involves a 2D interpolation and the re- 
sulting changes are discussed in the following. As before. 



we on ly consider data values fulfilling the criteria outlined in 
§[Z3l 

The binning causes a decrease of |0^| by ~3% (see Ta- 
ble [TJ. It also causes an increase of Bh ^sF^n ^y ~6%. The 
force-freeness remains basically the same, as do the relative 
occurrence of the vertical/horizontal field (Figure (T^/b) and 
the respective power distributions (Figure |7j:/d) remain simi- 
lar. Now it is also evident that the binning is only to a minor 
degree responsible for the diff^erenc es between t h e SPb Jn and 
HMI data. This agrees with, e.g., IWang et al.l (120091) who 
found that the scaling the data to a lower resolution does not 
significantly alter the results of their particular analysis. 

The HMI data hosts ~ 53% of l^^lsp^^g and ~ 86% of ^ sPong 
(see Table [T]). Naively, one would suspect that the diff'er- 
ence in |0^| arises from the diff'erent resolution limits of the 
two instruments, i. e., that the magnetic field is partially on 
scales which HMI cannot resolve. In context with the relative 
occurrence of |0^| discussed in Thalmann et al. ( 2012) , how- 
ever, this cannot be concluded since it was found that HMI 
vertical fields are on overall weaker than those of VSM data 
(which with a plate-scale of ~ 1 arc-second has a lower res- 
olution than HMI). Therefore, the diff'erent occurrence rates 
must have reasons besides the diff'erent resolution limit of the 
instruments. 

Preprocessing the SPorig data leads to a decrease of l^^lsPong 
by ~ 3% and and increase of Bh, sf^^ by ~ 40%, which is com- 
parable to the changes of the SPbin data due to preprocess- 
ing. The modification to the SPorig lower boundary data dur- 
ing solving for the force- and divergence-free field, as listed 
in Table [U is that l^^lsPon decreases by ~ 3% and Bh^sFon ^^~ 
creases by < 1%. Summarizing, we find a similar behavior of 
the modifications to the SPorig data during the force-free mod- 
eling and very similar values as we found for the SPbin data 
(compare the values listed in Table [TJ. 

The potential energy of the SPbin model is ~ 1 % lower and 
the total energy is ~ 2% higher than that of the SPorig model 
(see Table [2]). Hence, when basing the analysis on the SPbin 
data, we find and enhancement/reduction of the absolute en- 
ergy estimates £^niff/^pot on the order of the statistical error 
of the energy estimates itself. The excess energy, AE, is en- 
hanced by ~ 10%, i. e., also on the order of the statistical error. 
The relative amount of AE, however, remains approximately 
the same. 

Repeating the analysis of the magnetic connectivity within 
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Fig. 8. — (a) Squashing degree, Q, at a height of ~ 5 Mm above the NLFF lower boundary as determined from the closed SPorig model field lines, shown in the 
range x > 55 Mm. Field lines originating from locations where ^ > 100 within the ROI (outlined rectangular area) and which have both footpoints located on the 
lower boundary are shown in panels (b) and (c). The magnetic configuration when viewed along the surface normal is shown in (c). The gray-scale background 
reflects the vertical magnetic field, B,, of the NLFF lower boundary. Black and white contours are drawn at [-i-50mT,-i-100mT] and [-50mT,-100mT], 
respectively. The view from solar south (along y) is shown in (c). 




80 100 

X [Mm] 

Fig. 9. — NLFF lower boundary of the SPorig model in the range 
(55Mm<x,_y). The vertical magnetic field is shown as gray-scale back- 
ground. Black and white contours are drawn at [-i-50mT,-i-100mT] and 
[-50mT,-100mT], respectively. The locations where field lines originate 
from regions of 2> 100 within the ROI (rectangular outline) re-enter the 
lower boundary are color-coded according to the vertical magnetic field there. 

tlie SPorig model, we find that in total 228 magnetic field lines 
originating from locations of 2 > 100 within the ROI in Fig- 
ure [8^ connect back to the lower boundary (see Figure |8j),c). 
Again, we assume that those represent thin flux tubes with 
a cross section of d^ ~ 220^ km^ and assume the vertical 
field at footpoint from which we started the field line cal- 
culation determines its flux. Summation over the all con- 
sidered thin flux tubes gives an estimated total shared flux 
of I0JI = 1.2 X 10^^ Mx which comprises about 1% of the un- 
signed vertical flux of the NLFF lower boundary of the SPorig 
model. This is almost identical to what was found for the 
model based on the SPbm data, as is the connectivity pat- 
tern (compare Figures|9]and|5])). Considering all field lines 
starting from a pixel location on the NLFF bottom boundary 
in the range (60 Mm < x, j) where B^>2crB.^ and Bh,>2(TBh 
or B^ > 2(Tfi, and Bh < ^o-Bh ^ we find an identical relationship 
of /pLand /zmax (Figure [T0| as was found for the SPbm model 
(compare Figure (Q)). The SPorig model field is found to be 
even more vertical (< B^IBu > ~ 0.97) and, again, the field 
lines carrying the strongest vertical fluxes are those with an 
intermediate length and apex height. 

4. DISCUSSION AND CONCLUSIONS 

Non-identical photospheric vector magnetic fields inferred 
from polarization measurements of diff'erent instruments used 
as input for nonlinear force-free (NLFF) coronal magnetic 
field models directly translate to substantial diff'erences in the 
model outcome. To quantify these, we performed force-free 
magnetic field modeling using active-region vector magnetic 
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Fig. 10. — Relationship of the length of closed model field lines, /fl, to their 
apex height (/zmax) for the SPorig model. Field lines are calculated starting 
from every pixel location in the range (60 Mm < x, y) where B^ > 2o-b, and 
Bh,> 2o-Bh ^^ where B^ > Ictb. and Bh < ^ctb,^ . The color-code reflects the 
value of the absolute vertical field there. 

field data based on measurements of polarization signals by 
the SDO/HMI and Hinode SOT/SP. The possible causes of 
the diff'erences of the data products itself, including the, e. g., 
diff'erent intrinsic nature and sensitivity of the instruments, the 
temporal evolution of the photospheric magnetic field during 
the ongoing scanning times or the inversion techniques used 
to infer the magnetic field vector from the measured polar- 
ization signals, were out of the scope of this work. Our aim 
was to compare force-free model results based on data from 
these two instruments so we applied all data preparation and 
modeling routines in exactly the same way to both data sets. 

Force-free coronal magnetic field modeling is computation- 
ally expensive and high-resolution data is sometimes binned 
to a lower resolution in order to shorten the computational 
time. We, therefore, binned the original-resolution SP (SPorig) 
data to the resolution of HMI which allowed us, besides (i) to 
quantify the deviations of the force-free modeling outcome 
due to the usage of data of different instruments (HMI and 
SP), also (ii) to investigate the eff'ect of binning the input data 
to a lower resolution on the model outcome, by subsequent 
comparison of the models based on the binned SP (SPbm) data 
and SPorig data. We did not intend to mimic the diff'erent spa- 
tial resolution of the instruments by binning of the vector data 
since, as pointed out by Leka & Barnes (2012) , any kind of 
binning does not account properly for resolution effects. In- 
stead, in practice, a binning of high-resolution data prior to the 
NLFF modeling is meant to result in feasible computational 
dimensions and to allow for near real-time modeling. 

We used HMI and SP vector maps of active region 11382 



Thalmann et al. 



on 201 1 December 22 and found considerably higher vertical 
magnetic flux and average horizontal field in the SP data. This 
diff'erence of the vertical flux was found to be much larger 
than the modifications to the data due to application of our 
force-free modeling routines. Also the modifications to the 
vertical flux due to binning of the SPong data to the resolution 
of HMI were found to be small compared to the inequality of 
detected flux by the two instruments. 

Unequal estimates of the magnitude and orientation of lon- 
gitudinal and transverse fields based on the inversion of circu- 
lar and linear polarization signals measured by the two diff'er- 
ent instruments, yield diff'ering estimates of the vertical and 
horizontal magnetic field in a local coordinate system. The 
unequal amount of vertical magnetic flux (the projected SP 
data hosting about two times the unsigned vertical flux of 
HMI) then directly translates to a considerable discrepancy 
in the absolute estimates of the energy content of the consid- 
ered coronal volume: models based on the SP data hold about 
twice as much energy as does the HMI model. By tracing 
magnetic field lines in the half- space above the model lower 
boundaries their connectivity was investigated. 

The SP and SPorig models revealed a closely matching pho- 
tospheric footprint of two neighboring coronal connectivity 
domains. A similar footprint was found in the HMI model but 
locally shifted by up to several Mm. The same applies to the 
locations where field lines re-enter the NLFF lower boundary : 
while the location of the footpoints of the SP and SPong mod- 
els coincide, those of the HMI model are displaced by up to 
ten Mm. Moreover, on overall, the SP and SPong model fields 
tend to be more vertical than the HMI model field which is a 
direct consequence of the relative vertical and horizontal field 
distribution given by the instrument data. 

However, the models also showed great similarities: rel- 
ative estimates like the fraction of energy in excess over a 
potential field (about 20% of the total energy content) or the 
fraction of vertical flux shared by two neighboring connectiv- 



ity domains (about 1 % of the total vertical flux of the active 
region) agree very well. Also the overall field-line geometry 
was found to be comparable: the length of all closed field lines 
in the model volumes was found to be more or less linearly re- 
lated to the apex height. Common to the model outcomes is 
also that the shortest and lowest as well as the longest and 
highest field lines carry least vertical flux. 

In conclusion, caution is needed when analyzing the coro- 
nal magnetic field and its connectivity with the help of force- 
free magnetic field models based on the vector magnetic field 
products of diff'erent instruments, made available to the com- 
munity. Relative estimates and the overall structure of the 
model magnetic fields might indeed be reliable while absolute 
estimates might only be so concerning their order of magni- 
tude. Moreover, binning of the magnetic vector data to a lower 
resolution prior to the force-free modeling results only in little 
diff'erences in the model outcome, small compared to the re- 
markable deviations when basing the modeling on data from 
the two diff'erent instruments, HMI and SP. 
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